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Abstract 

The design of heat exchangers traditionally focuses on the known con- 
straints of the problem such as inlet and outlet temperatures, flow rates, and 
pressure drops. This leads mainly to a sizing problem where the designer 
must select surfaces, flow configuration, and materials to meet the mini- 
mum design objectives. An alternate approach based on an acceptable level 
of thermodynamic irreversibility (entropy generation) has been proposed. 
When the entropy generation level has been set, the geometric parameters 
of the heat exchanger can be determined. The design of a plate-fin type, 
gas-to-gas recuperator for a regenerative open Brayton cycle has been used 
as a demonstrative device. The resulting heat exchanger designs are then 
examined to determine what caused the differences and why either method 
should be preferred over the other. 
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Chapter 1 

INTRODUCTION 



1.1 Background 

Energy conservation is a topic that has recently received considerable atten- 
tion. Available work has been recognized as a valuable commodity and its 
destruction in seemingly otherwise efficient engineering processes has also 
seen a renewed increase in interest. The use of second law analysis and ther- 
modynamic irreversibility minimization have been proposed as techniques 
that should become an integral part of the design of engineering processes 
and components [l]. It is not apparent what the relationship of traditional 
methods of design and these irreversibility minimization techniques will be. 

1.2 Conventional Design Procedures 

The design of a component, such as a heat exchanger, generally involves the 
sizing of that component to meet specified performance parameters within 
known constraints. The heat exchanger area and volume are usually de- 
signed to be the minimum required to meet the specifications as this will 
also usually be the best design from an ecomonic aspect. The fluid inlet and 
outlet temperatures, flow rates, and pressure drops are usually specified and 
it is the designer’s task to determine construction type, flow arrangement, 
materials and surfaces to meet those requirements [2]. The effect of the 
component design on the overall system performance should also be consid- 
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ered to ensure the tradeoffs performed within its design are still valid when 
integrated into the system. 

1.3 Irreversibility Minimization 

Any heat transfer process is generally accompanied by thermodynamic ir- 
reversibility or entropy generation. The entropy generated is in direct pro- 
portion to the amount of useful work dissipated in that process. If this 
irreversibility can be minimized within a particular component of a power 
cycle, the useful power output of the cycle should increase [3]. Bejan has 
shown how the entropy generation rate can be reduced in a counterflow gas- 
to-gas heat exchanger [4], and it is this method that will be used to create 
designs that will be compared to conventionally designed heat exchangers, 
and their impact on overall cycle performance. 
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Chapter 2 

SELECTION OF 
BRAYTON CYCLE 
PARAMETERS 



2.1 Cycle Components 



A regenerative open Brayton cycle was chosen as the vehicle to test the 
different design methods for the recuperator. A schematic of the cycle com- 
ponents is shown in Figure 2.1. 




Regenerative Cycle Schematic 

Figure 2.1 
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A compressor polytropic efficiency was assumed and used throughout 
all calculations, t] cp = 0.95. Similarly, a turbine polytropic efficiency was 
assumed, rj ct = 0.90. For ease of calculations, only one fluid was used, 
i.e., air was modeled as a perfect gas with c p = 1.0 kj/kg — ° K and R = 
0.287 kj/kg— ° K both assumed constant. The combustor, then, was treated 
as a perfect heat transfer device, but was not an injection point for fuel. 



2.2 Cycle Parameters 

A reasonable compressor inlet temperature was selected, T\ = 300° if, along 
with an inlet pressure of Pi = 1.126 X 10 5 N/m 2 . A nominal value of turbine 
inlet temperature was also selected, T 4 = 1300° i£\ The selection of an 
appropriate pressure ratio required the following consideration. 

The addition of a recuperator to a simple open Brayton cycle causes a 
shift in the pressure ratio for maximum cycle efficiency from that of the non- 
regenerative cycle. A heat exchanger effectiveness of e = 0.8 was selected 
and cycle efficiency without losses was calculated for pressure ratios varying 
from 4 to 10. A sample calculation is contained in Appendix A, results are 
plotted in Figure 2.2. 



Cycle Efficiency vs Pressure Ratio 
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Based on these results, a pressure ratio of PijPi = 5 was selected. With 
this pressure ratio determined, the maximum cycle efficiency without losses, 
7700 , could be calculated. This quantity will be used for comparison purposes 
and was computed as r}^ — 0.5760. 



2.3 Heat Exchanger Characteristics 

A plate-fin heat exchanger surface was selected from those cataloged in ref- 
erence [5]. Its characteristics and dimensions are shown in Figure 2.3. It has 
an uncomplicated geometry and a reasonably wide range of Reynolds num- 
bers in which consistent results can be expected. The material is aluminum, 
with a thermal conductivity oikj = 190 w/m — ° K . 




Fin pitch = 2.0 per in = 78.74 per m 
Plate spacing, b - 0.750 in = 19.05 x 10" 3 m 
Fin length = 12.0 in = 304.8 x 10° m 

Flow passage hydraulic diameter, 4 r h = 0.0474 ft = 14.453 x 10' 3 m 

Fin metal thickness* 0.032 in, aluminum = 0.813 x I0' 3 m 

Total heat transfer area/volume between plates, 0 = 76.1 ft 2 /ft 3 = 249.672 m 2 /m 3 

Fin area/total area = 0.606 



Plain Plate-Fin Surface 2.0 

Figure 2.3 
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Chapter 3 

CONVENTIONAL 
DESIGN METHOD 



3.1 Introduction 

This chapter outlines the conventional design procedure that was used to 
size the single-pass, counterflow gas-to-gas heat exchanger. This procedure 
was performed for various combinations of heat exchanger effectiveness and 
specified pressure drops. A complete numerical example is provided in Ap- 
pendix B. Final results for this conventional method are shown in Table 
3 . 1 . 

3.2 Conventional Methodology 

For a single-pass counterflow design, when the core dimensions on one side 
are fixed, the dimensions on the other side are also known. This means the 
design is driven by the side that has the more stringent A P/P requirement. 
The method of determining the controlling side is contained in Appendix 
B. In every case considered, the hot side of the heat exchanger was the 
controlling side and so its requirements drove the design. 

The following is a step-by-step heat exchanger sizing design procedure 
which closely follows that outlined in reference [2]. 

1. For the given heat exchanger effectiveness, determine the fluid outlet 
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temperatures. Calculate the fluid mean temperature on each side and 
evaluate fluid physical properties p 0 , (l/p) m , and 

2. Determine N tu for the exchanger and then N tu for each side. The in- 
fluence of longitudinal heat conduction is ignored in this first iteration 
of design. 

3. Estimate hot and cold side pressure drops, select an appropriate Nr 6 , 
and then a value of jj f from Figure 2.3. 

4. Calculate mass velocity, G, from information in steps 1-3 and the cor- 
responding value of A P/P. 

5. Calculate Nr c and determine values of j and f from Figure 2.3. 

6. Compute heat transfer coefficient, h; temperature effectiveness of the 
fins, rjf\ and the total surface temperature effectiveness, rj Q . 

7. Calculate heat transfer area, A\ minimum free flow area, A 0 ; heat ex- 
changer frontal area, Af r \ flow length, L; and heat exchanger volume, 
Vol. 

8. Compute A P/P from known conditions and calculated parameters. 

3.3 Analysis 

For each case show in Table 3.1, cycle efficiency with losses, 77, was calculated 
and the cycle efficiency degrade, — 77, is listed for each case. A correlation 
between heat exchanger effectiveness, e, total allowed pressure drop, k, cycle 
efficiency degrade, and heat exchanger volume was sought. 

The relationship between cycle efficiency degrade and effectiveness at 
various levels of k is shown in Figure 3.1. The increased level of degrade 
at lower values of effectiveness is as expected, and for a particular value of 
effectiveness, the degrade increases with increasing pressure drops. 

The effects of varying effectiveness and pressure drops on heat exchanger 
volume is depicted in Figure 3.2. The higher the effectiveness, the longer 
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and therefore, larger the heat exchanger. Similarly, the greater pressure 
drop for a particular effectiveness also increases the volume. 

Figure 3.3 was derived by extracting values of effectiveness at nominal 
values of r ] OQ — 77 and plotting them on the appropriate k curve of Figure 3.2. 
This resulted in curves of constant — r\ that have a minimum volume 
at a particular value of k and e. It is emphasized that these “minimum” 
volumes were obtained at constant values of cycle degrade with varying 
pressure drops and effectiveness. These “minimum” volumes are plotted 
against effectiveness in Figure 3.4 with the respective k values annotated on 
the graph. 
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Case # 


e 


k 


43 

8 

1 

43 


L(m) 


Vol (m 3 ) 


C-l 


1-0 


0.02 


0.0043 


— 


— 


C-2 


1.0 


0.04 


0.0089 


— 


— 


C-3 


1.0 


0.06 


0.0137 


— 


— 


C-4 


1.0 


0.08 


0.0188 


— 


— 


C-5 


0.9 


0.0 


0.0449 


— 


— 


C-6 


0.9 


0.0240 


0.0498 


13.69 


3.69 


C-7 


0.9 


0.0388 


0.0550 


13.82 


2.64 


C-8 


0.9 


0.0558 


0.0603 


14.04 


2.19 


C-9 


0.9 


0.0739 


0.0660 


14.65 


1.98 


C-10 


0.8 


0.0 


0.0833 


— 


— 


C-ll 


0.8 


0.0192 


0.0886 


6.15 


1.12 


C-12 


0.8 


0.0367 


0.0941 


6.56 


0.85 


C-13 


0.8 


0.0550 


0.0999 


6.72 


0.71 


C-14 


0.8 


0.0740 


0.1060 


7.15 


0.65 


C-15 


0.7 


0.0 


0.1165 


— 


— 


C-16 


0.7 


0.0188 


0.1221 


3.75 


0.53 


C-17 


0.7 


0.0365 


0.1278 


3.95 


0.39 


C-18 


0.7 


0.0565 


0.1338 


4.17 


0.34 


C-19 


0.7 


0.0734 


0.1401 


4.32 


0.31 


C-20 


0.6 


0.0 


0.1455 


— 


— 


C-21 


0.6 


0.0186 


0.1512 


2.52 


0.29 


C-22 


0.6 


0.0376 


0.1571 


2.69 


0.22 


C-23 


0.6 


0.0557 


0.1632 


2.85 


0.19 


C-2 4 


0.6 


0.0764 


0.1695 


2.98 


0.17 



Table 3.1: Conventional Method Results 
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Efficiency Degrade vs Effectiveness 

Figure 3.1 
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/ olume vs Effectiveness 

Figure 3.2 
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Volumes at Constant Cycle Degrade 

Figure 3.3 
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Minimum' Volume vs Effectiveness 

Figure 3.4 
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Chapter 4 

MINIMUM ENTROPY 
DESIGN METHOD 



4.1 Introduction 

The use of irreversibility concepts in the design of the gas-to-gas counterflow 
heat exchanger will be presented in this chapter. For the same mass veloci- 
ties that were used in Chapter 3, the heat exchanger geometric parameters 
can be determined to arrive at minimum levels of entropy generation. The 
outputs of this design procedure are heat exchanger effectiveness and the 
pressure loss suffered on each side. For a more detailed explanation of this 
procedure, the reader is directed to references [4] and [6]. 

4.2 Background 

Since the method deals with the entropy generation within the recuperator, 
it is appropriate to begin with an entropy flux analysis on the control volume 
of the recuperator in Figure 4.1. The heat transfer from the outer walls is 
assumed to be negligible. The entropy generation rate can then be written 
as 

S :=z $tn)l ~b T7l2 ^ tn )2 (4.1) 
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Stream 1 
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Stream 2 



t 2 




Recuperator' Control Volume 

Figure 4.1 



Substituting expressions for the ideal gas entropy changes based on pressure 
and temperature relationships, this becomes 

s = C min [ln{T l0 ut/T{) + {R/CjMPi/Piout)] 

+C max {ln{T 2otlt /T 2 ) + (R/C p ) 2 ln(P 2 /P 2out )} (4.2) 

The number of entropy production units, N 3 , is defined as 

N, = S/C max (4.3) 

Writing expressions for the first law of thermodynamics of the control vol- 
ume, 

c min (Tj — Tw) + C m i n (T 2 — T 2 out) = 0 (4.4) 

and heat exchanger effectiveness, 

€ = {Tiaut ~ Tx)/{T 2 - Tx) (4.5) 



24 



\ 



permits the elimination of Ti out and T 2ou t in equation (4.2) which allows the 
rate of entropy production formula to be nondimensionalized as 




A special form of this last expression can be obtained in the case of nearly 
ideal heat exchangers. When the stream-to-stream AT’s and frictional AP’s 
are held small, two inequalities can be applied. 

• l-e << 1 (4.7) 



(A P/P)i,2 « 1 



Applying these and expressing e in terms of Ntu yields 



N s 



C m in 

C m ax 



In 
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Cmin 

Cmax 




(4.8) 



( Cmin Y(i 
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! i 1 ~ ft. 


) ! 


exp [-iV fu (l - 


Ts 

O O 


\c max ) \ 


Cmax ) 


i - (l 


-ft) 


1 - £Wezp I -N tu 

'-'max L 


(l - pain.) 

\ '-'maz / . 




(4.9) 



The first two terms can be characterized as a contribution from capacity 
rate imbalance, the third represents contribution due to finite N tu , and the 
last two show fluid friction effects. 

Requiring the streams to be balanced and applying the calculus of limits 
as C m in — > C max , equation (4.9) can be reduced to 



N s 



Cmin j T 2 

In — 

Cmax 3 \ 



+ 



In 




C m in 

Cmax 



Ti 

T 2 



)] 
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i 
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C , (410 > 



When the overall N tu is expressed in terms of N tu for each side of the heat 
transfer surface, 

1 — 1 _j_ Cmin 1 

iV* u iV^ui 

Cmax Ntu 2 

the number of entropy generation units can be divided into three contribu- 
tions 

N 3 — N 3 imbalance “f" N 3 \ + N 3 2 (4.12) 



Applying the limit C miri 
can be written as 

N 3 imbalance 




Cmax , the imbalance terms of equation (4.10) 



- 1 




C m in 

Cmax 




(4.13) 



which will clearly vanish when C m i n = C max . Since this is the case for the 
design here, the number of entropy production units for each side can be 
determined from 




(4.14) 



(4.15) 



4.3 Design Considerations 

The number of entropy generation units that have been deemed acceptable 
for a particular design is a function of the hydraulic radius, minimum free 
flow area and flow length. It is necessary, then, to express N s in terms of 
these and other known flow parameters. Using the definitions of Nt u and 
friction factor for each side 

Ntui, 2 = ( L/r h )N st (4.16) 



(A P/P)!, 2 = f(L/r h )G 2 /(2pP) (4.17) 
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(4.18) 



in equations (4.14) and (4.15) leads to 

CLT 

where 



CL i = 1 



^2 — Cmin/ G n 



b\ — Cmin/Cv 



b 2 = 1 

Bi = (R/C p )i 

B 2 = (R/C p ) 2 




g = G/(2pP)* 



(4.19) 



(4.20) 



(4.21) 

(4.22) 



(4.23) 



As can be seen in equation (4.18), the ratio L/r^ performs a trade-off func- 
tion, i.e., for a fixed g and Nr 6 , there will be an optimum L/r^ which results 
in a minimum N $ i t2 . When the function N s is minimized, the ( L/r h ) opt 
is given by 



{L/ r h)opt 



2 / r/N at 

G \Bf(l/p c P c + l/p h P h ) 



(4.24) 



4.4 Minimum Entropy Design Methodology 

The method used is based on optimizing (L/r^). Since is set by the selec- 
tion of the plate-fin surface, this yields a length that will produce minimum 
entropy for the given hydraulic radius. 
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The following outlines the methodology employed. A complete numerical 
example is shown in Appendix C. Final results for all cases are tabulated in 
Table 4.1. 

1. For the cases listed in Table 3.1 that resulted in a heat exchanger 
volume, use calculated G and Nr € to determine ( L/rh) opt . 

2. Compute resultant heat exchanger effectiveness and hot and cold side 
pressure drops. 

3. Calculate N Sfi , N Sc , and N Stotal , 

4. Calculate the heat exchanger volume for each case. 

5. Determine the cycle efficiency degrade, — rj. 

4.5 Analysis 

Figure 4.2 shows the nearly linear relationship of cycle efficiency degrade and 
heat exchanger effectiveness. There is not the clear definition of different 
values of pressure drop as was observed for the conventional method in 
Figure 3.1. In fact, there is even a wider range of k values here, from 
k = 0.056 for case E-6 to k = 0.269 for case E-24, and yet the results are 
almost linear across that range. 

The volume plotted in Figure 4.3 is that volume obtained by minimizing 
entropy generation in the heat exchanger and should not be interpreted as 
a minimum volume, but rather as an “optimum” volume for minimizing 
entropy generation. 

To achieve a better basis for comparison of the two methods, the entropy 
generated by the various designs must be discussed. 
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Case # 


L opt {m) 


Vol(m 3 ) 


k 


€ 


Voo-V 


n 3 


E-6 


32.03 


8.63 


0.0563 


0.955 


0.0347 


0.0325 


E-7 


25.86 


4.93 


0.0727 


0.945 


0.0440 


0.0417 


E-8 


22.04 


3.44 ' 


0.0876 


0.935 


0.0528 


0.0503 


E-9 


20.12 


2.72 


0.1014 


0.926 


0.0606 


0.0583 


E-ll 


24.35 


4.46 


0.07658 


0.942 


0.0461 


0.0441 


E-12 


19.08 


2.46 


0.1069 


0.922 


0.0641 


0.0614 


E-13 


16.11 


1.07 


0.1320 


0.908 


0.0784 


0.0760 . 


E-14 


15.08 


1.38 


0.1559 


0.896 


0.0917 


0.0897 


E-16 


19.80 


2.77 


0.996 


0.927 


0.0596 


0.0573 


E-17 


15.40 


1.53 


0.1422 


0.902 


0.0844 


0.0823 


E-18 


13.19 


1.07 


0.1778 


0.884 


0.105 


0.1028 


E-19 


12.11 


0.86 


0.2058 


0.872 


0.1211 


0.1189 


E-21 


16.98 


1.93 


0.1252 


0.913 


0.0741 


0.0718 


E-22 


12.87 


1.04 


0.1799 


0.882 


0.1067 


0.1048 


E-23 


11.52 


0.76 


0.2251 


0.862 


0.1335 


0.1309 


E-24 


10.47 


0.60 


0.2685 


0.847 


0.1596 


0.1541 



Table 4.1: Entropy Generation Method Results 
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Eff. Degrade vs Min Entropy Effect. 

Figure 4.2 
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iin Entropy Volume vs Effectiveness 

Figure 4.3 
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Chapter 5 

COMPARISONS AND 
RESULTS 



5.1 Introduction 

To this point, the minimum entropy method results have been compared to 
the conventional method results in terms of conventional method standards, 
i.e., pressure drops, efficiencies, and effectiveness. It is necessary to compare 
these methods in terms of the minimum entropy method, i.e., the entropy 
generated in each design. 



5.2 Conventional Method Entropy Generation 

The number of entropy generation units for each design can be computed 
from previously determined information and equations (4.14) and (4.15). 
An example of these calculations is contained in Appendix D. Results for 
the conventional cases are listed in Table 5.1. The large deviation in N s and 
Voo ~ *7 for cases C-21 through C-24 are a result of violating the assumption 
that 1 — e << 1. These conventional cases will be discarded from further 
consideration. However, the minimum entropy cases at corresponding mass 
velocities are still valid. 

When the entropy generation units are plotted against the degrade in 
cycle efficiency, Figure 5.1, the dependence on pressure drop is still present. 
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The grouping of data points is caused by the discrete changes in heat ex- 
changer effectiveness, decreasing from left to right. 

Figure 5.2 was created by using the same values of ry ^ — ry that were used 
for Figure 3.4. The value of N s at each of these cycle efficiency degrades 
was read off at a particular value of k. Plotting these “minimum” volumes 
against entropy generation rate exhibits the general trend that the larger 
the volume of the heat exchanger, the lower the total entropy generated. 
This is in consonance with the results predicted by reference [4]. 

5.3 Minimum Entropy Generation 

The cycle efficiency degrade vs. entropy generated for the minimum entropy 
method is similar in overall shape to that of the conventional method, but, 
as Figure 5.3 shows, the relationship is nearly linear and does not show the 
discrete pressure drop differences. 

Figure 5.4 shows that the “optimum” volumes determined by (L/r^opt 
follow essentially the same path as that of the conventional method. For a 
direct comparison, Figures 5.2 and 5.4 are plotted together on Figure 5.5. 

To ensure that the (L/rh) op t expression, (4.24), was giving minimum 
entropy results, a series of calculations was performed at ( Ljr values on 
both sides of the optimum for case E-9. Figure 5.6 shows that the minimum 
entropy solution was, in fact, being determined. It is also noted that the 
magnitudes of N s and ry^ — ry are essentially the same over this range of 
heat exchanger effectiveness. 

5.4 Comparisons 

The calculations performed on case E-9 to prove minimization also yielded 
heat exchanger volumes at a constant mass velocity G. These additional 
data points are plotted on Figure 5.7. The dash-dot line connecting them 
is a locus of recuperator volumes along a constant G path. This shows why 
the volumes at minimum N s or ry^ — ry along a constant G path are larger 
that the miminum volume of a given N 3 or ry - ry. 
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Figure 5.8 shows a comparison between conventional “minimum” vol- 
umes and minimum entropy “optimum” volumes. The minimum entropy 
method yields volumes that are larger than the conventional method for the 
same cycle efficiency degrade. 
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Case # 


C3 

8 

1 

cs 


N, 


C-l 


0.0043 


— 


C-2 


0.0089 


— 


C-3 


0.0137 


— 


C-4 


0.0188 


— 


C-5 


0.0449 


— 


C-6 


0.0498 


0.0445 


C-7 


0.0550 


0.0494 


C-8 


0.0603 


0.0550 


C-9 


0.0660 


0.0610 


C-10 


0.0833 


— 


C-ll 


0.0886 


0.0900 


C-12 


0.0941 - 


0.0967 


C-13 


0.0999 


0.1036 


C-14 


0.1060 


0.1108 


C-15 


0.1165 


— 


C-16 


0.1221 


0.1505 


C-17 


0.1278 


0.1584 


C-18 


0.1338 


0.1669 


C-19 


0.1401 


0.1748 


C-20 


0.1455 


— 


C-21 


0.1512 


0.2308 


C-22 


0.1571 


0.2405 


C-23 


0.1632 


0.2501 


C-24 


0.1695 


0.2607 



Table 5.1: Conventional Method Entropy Generation 



35 



Efficiency Degrade vs Entropy Generated 

Figure 5.1 
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Minimum' Volume vs Entropy Generated 

Figure 5.2 
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Figure 5.3 
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Optimum' Volume vs Min Entropy Gen. 

Figure 5.4 
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Volume vs Entropy Comparison 

Figure 5.5 
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Entropy Generation Unite 

minimum entropy + conventional method 



Case E— -9 (L/V h ) Optimization 

Figure 5.6 
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Figure 5.7 
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Figure 5.B 
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Cycle Efficiency Degrade 

minimum entropy ■+ conventional method 



Chapter 6 

CONCLUSIONS 



6.1 Discussion 

The nature of the results and conclusions drawn for the example cycle would 
be the same for other cycles and recuperator surfaces. For high effective- 
ness recuperators, the counterflow arrangement is dictated. Here with the 
same finned surface on both sides, the (AP/P) c /(AP/P)/ l ~ (P h /P c ) 2 or 
approximately 1 / 25 . If a closer plate spacing were used on the clean air cold 
side more A P/P would appear on the cold side for any total S P/P. This 
would result in even smaller minimum volumes at any given N s or — r\. 

6.2 Conclusions 

In selecting a gas turbine recuperator, essentially the same results are ob- 
tained either by the conventional method of minimizing degradation of cycle 
efficiency, 77 oo - 77, or by minimizing entropy generation, N s . 

The selection of an “optimum” recuperator requires that care be exer- 
cised in the selection of the path along which the optimization is calculated. 
It has been shown that the volume of the recuperator for minimum N s or 
77 00 “ 7 ] along a constant mass velocity path is much greater than the true 
minimum volume at any 77 oo — 77 or N s . 
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Appendix A 

Sample Calculation for 
Effect of Recuperator on 
Cycle Efficiency 



e = 0.8, r = - 



Wnet 

»7co = T 

Q 



W t + iu c 

Q 



u t — m(/i 4 — /15) = mc p (T 4 — T5) 



cj c = rn{h\ — /12) = mc p (Ti — T{) 



Q = m(/i 4 - /i 3 ) = mCp(T 4 - T 3 ) 

Turbine: 

74 _ 

57 " VA7 



t 5 




(^)(^) 



t 4 



-1 
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Compressor: 



Recuperator: 



T s = 



( 6 )(' 9 )(' 



-1 



1300 



= 820.1°fc 



T 2 / P 2 \ (’ipcH-i) 

Ti = \TJ 

n = r,^)^ 

T 2 = (300) ( 6)^07 =51 4.2° k 



T2-T2 
6 T 5 - T 2 

T 3 = (e)(T 5 -T 2 )+T 2 



T 3 = (.8) (820.1 — 514.2) + 514.2 = 695.2° fc 
Cycle Efficiency: 

(T 4 - T s ) + (Ti - T 2 ) _ (1300 - 820.1) + (300 - 514.2) 
1100 ~ T 4 - T 3 “ 1300 - 695.2 



0.4651 
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Appendix B 

Conventional Design 
Method Numerical Example 



Before calculations can be performed, two points need to be clarified: 1) the 
method of determining how loss of turbine work due to pressure drops in 
the system is distributed and; 2) the method for determining which side of 
the heat exchanger has the more stringent A P/P requirement. 

Turbine work is given by the expression 

w t = mr) t c p T h 




w t 



rnr)tC p T h 








’i+^' 


\Pk) 


L 1 "^] 



Wt 



mr) t c p T h 



Wf 

rnrjtCpTh 



last term, 








~ 1 - 1 


'PS 


s 


i 1+ pS 


- 1 1 


\Ph> 


1 


- I" 1 


'PS 


r 


r tJ _ m Pk 


\Ph; 


/ 



+ 7T' + 



so that loss of turbine work can be expressed as 

mr, t c p T h 1 \Pj V P h Pc J 
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and it is clear from Figure B.l that determines the loss of turbine 

work irregardless of how it is distributed. Introducing the definitions r = Q 




specific entropy 

Temperature - Entropy Diagram 
Figure B.l 



and 9 — (SM , it can be shown that for a counterflow exchanger where 

V ■‘c / avg 

A e — A h and L c = Lh, 

APh/Ph _ 9 

APc/Pc r 2 

Using a slightly different form of the last turbine work expression, 



Aw t 


1 


\ APh/Ph 


+ 1 


AP C 


mritCpTh 




[ap c /p c 


Pc 


Aw t 


_ * i 


(? +1 ) 


A P c 




rmitCpTh 


7 r 7 


Pc 
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Ati>t 
ti ^ 



1 ( 0/r 2 + 1) A P c 



Awt _ {e/r 2 + l)AP c 

w t ~ 7 r^-l P c 

Introducing = k, expressions for hot and cold side pressure 

drops are, 



A Pk = A:(0/r 2 ) 

Ph 1 + 0/r 2 



AP C _ A: 

Pc' ~ 1 + 0/r 2 



The numbered steps that follow correspond to those outlined in Section 
3.2. The input design parameters used for the example calculations were 
e — 0.8 and k — 0.04. 



B.l Temperatures and Fluid Properties 



T x = 300° k 



so that 



Ti 

Ti 




npci 



r 2 



t i 




npc“t 



T 2 = 300(5) = 486.8° A: 



Pa 

Ps 



Pi 
P i 



(1 - A:) = 5(1 - .04) = 4.8 
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Ti = ( 

Ts \ 



*) 

Pj 



Hpt 






Ts = 



(%r 






-i 



T 4 



Ts = 



'(4.8p)(^) 



-l 



1300 



= 868.5° k 



so that 



Similarly 



e = 



T3-T2 

Ts-T 2 



T 3 = e(T 5 -T 2 ) + T 2 



Ts = (.8)(868.5 - 486.8) = 792.2°*: 



T 6 = Ts-e{Ts-T 2 ) 



Tq = 868.5 - (,8)(868.5 - 486.8) = 563.1°*: 

Since this is a counterflow design with C* > 0.5, the best choice for 
average temperatures is the arithmetic average on each side 

T 5 + T 6 868.5 + 563.1 



'^'hrn — 



= 715.8°*; 



JU = = 639.5’* 

Inlet, outlet and average densities are next computed. Pv = RT or 
v = RT j P but v = 1/p, so p — P / RT 

1.126x10 s N/m 2 , 3 

* = WWs ~ *0(868.5°*) = 
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P 6 



1.126 x 10 5 
(287)(563.l) 



0.6967 kg/m 3 




\ (— + — ) = 1.824m 3 /^ 

2 V P5 P6 / 



P2 



5.63 x 10 5 
(287)(486.8) 



4.030A:<7/m 3 



P 3 



5.63 X 10 5 
(287)(792.2) 



2A76kg/m 3 




0.3260m 3 /kg 



B.2 iV (u 



For a couterflow heat exchanger, 

e _ 1 - exp[-jV tu (l - c*)] 

1 - c*eip[-iV tu (l - c*)] 

but for the balanced flow case (c* = 1) under design, 

Ntu, 

1 + 



or 



N t u = 



0.8 



= 4 



1 - e 1-0.8 

Since gas is the working fluid on both sides, it was estimated that both sides 
would have approximately the same surface resistance so that 



Ntu,» — Ntu t c — 2i\ft u — 8 



B.3 Estimated Pressure Drops 



Recall 




715.8 

639.5 



1.1193 
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and 



9 _ 1.1193 

(V5) 5 ' 



= 27.983 



= k{0(r>)_ = (0.04) (27.983) = 

P h 1 + 9/r 2 1 + 27.983 



A Pc 



0.04 



1 + 0 Jr 2 1 + 27.983 



= 0.00138 



From Figure 2.3, an estimate of j / f = 0.0035/0.0077 was obtained. 

B.4 Mass Velocities 

The mass velocity on each side is estimated from the relationship 
G = 



( 2g c Pt] 0 \ 


/AP/P> 


) - 


- V ( V P)mNp r J 


V A^tu > 


side / 



As a first approximation, r? 0 was assumed to be 0.90. 

" (2)(1.126x 10 s N/m 2 ) (0.9) \ / 0.03862\ /0.0035\ 



Gh — 

(1.824m 3 /kg)(0.1)* 

Gh = 17.586A:p/m 2 — s 

(2)(5.63 X 10 5 



!\ ( 0.0035 \ 
’/ V 0.0077/ 



Gr — 



G c — 



) (0.9) 'N / 0.00138 \ /0. 

if ) V 8 ) Va 



0035 \ 
0077/ 



NRe.h 
NRe,t i 



. (.3260) (0.7) J 

17.583 kg/m 2 - s 

G h D h (l7.583Jfcg/m 2 - s)(0.0474/t)(.3048m//t) 



M/i 

7481 



3.396 X 10 5 kg/m — s 
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_ (17 583) (.0474) (.304 8) = 

Re,c 3.159 x 10- 5 

Entering Figure 2.3 with these values yields 



f h = 0.0078 
j h = 0.00355 
f c = 0.0077 
j c = 0.0035 



B.5 Heat Transfer Coefficients and Fin Effective- 
nesses 

With the Colburn factor known, the heat transfer coefficient can be deter- 
mined from 

h = jGc p Pr * 

This result is used to compute the fin characteristic length, 



m = 



2h 

k f S 



The temperature effectiveness of the fin is next calculated from 

tanhml 



Vf = 



ml 



A good approximation for / is given by half the plate spacing minus the fin 
thickness. 

The total surface temperature effectiveness can then be calculated 

77o = 1 - (1 - r] f )Afl A 

where Af/A is a characteristic of the surface selected. 

(0.00355)(l7.586A:y/m 2 — s)(l.0kj/kg — ° k) 



h h = 



(0.7) 
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h h = 1 


f 9A9w/m 2 —° k 


m h - 


(2)(79.19 w/m 2 -°k) 1* 

( 190w/m — ° fc)(0.032m)(0.0254m/m) 


m h = c 


12. 024m -1 


' = 1 


( 0.750 . . \ 

— - — in — ,032m J 0.0254m/m 


l = f 


1.71 x 10 -3 m 


nf.h = ■ 


!an/i(32.024)(0.71 X 10 -3 ) 
(32.024)(0.71 X 10- 3 ) 


nf,h = ( 


).9749 


no ,h = i 


L- (1 -0.9749) (.606) 


no,h = ( 


).9848 


K = • 


C0.00 3 55)(17.5S 3 )(1.0) _ _ 

(0,7)1 


m c — 


(2)(78.06) 1* ? , 

L(190)(0.032)(0.0254)J 


n/,c = ■ 


tanh(31.79)(8.71 X 10- 3 ) _ 
(31.79)(8.71 X 10- 3 ) 


*70,0 — 


L - (1 - ,9752)(.606) = 0.9850 



The overall heat transfer coefficient can now be determined neglecting foul- 
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ing resistances and wall thermal resistances, 



1_ 

U 

1 

u 



1 1 

+ 



(r]oh)h ( f /o^i)c 



+ 



(0.9848)(79.19) (.9850) (78.06) 



= 0.02583 



U = 38.714 w/m 2 -° k 



B.6 Dimensions 

The heat transfer surface area is computed from 
N tu C 



A = 



U 



_ (4)(3600kg/hr)(1.0kj‘/kg -° k) _ 2 

38.714u;/m 2 —°k 

Minimum free flow area is given by 

= (W/O) = = 0.05686m ! 

v ' ’ 17.586^/m 2 - s 

In order to determine frontal area, a must be computed from the relationship 

dfiD h /A ■ 
d d* d + It 

Here a parting plate thickness of 0.4mm was assumed, so 

(0.75m)(76.1/t 2 //*3)(0.0474/4/t) 



a — 



0.75m + 0.75 in + 2(.4mm)(.l/2.54 m/mm) 



a = 0.44162 



Af r = Ao/c 



0.05686m 2 , 

At r = = 0.12876m 2 

;r 0.44162 
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The flow length is given by 

(D h A\ 



L = 



L = 



V4A 0 ) 

(0.0474/t) (0.3048m/ /t)( 103. 321m 2 ) 
(4) (0.05686m 2 ) 



L = 6.56m 

The total volume of the heat exchanger is calculated by 
Vol =(A fr )(L) 



Vol 



(0.12876m 2 ) (6.56m) = 0.845m 3 



B.T Calculated Pressure Drops 

For this level of design, entrance and exit effects are not considered, so only 
core friction effects will be used. 

A P 

= f(L/r h )(l/p) m G>/2P 



APh (0.0078) (6.56) (! .824) (!7.586) 2 _ 

P h (2) (0.0474/4) (.3048) (5.63 x 10 5 ) 



Similarly 



A P c 



= 0.00125 



so that 



AP 

E— = 0.03674 



The cycle efficiency with these input parameters is 

{t 4 - r 5 ) + (T\ - r 2 ) 
v t a -t 3 



(1300 - 868.5) + (300 - 486.8) 
n ~ (1300 - 792.2) 
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1 



0.4819 



and the efficiency degrade is given by rj <*, — r). 

0.5760 - 0.4819 = 0.0941 
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Appendix C 

Minimum Entropy Design 
Method Numerical Example 



This example uses the data from case C-12 of Table 3.1, 
G = 12.456 kg/m 2 - s, j = 0.00355, and / = 0.0078. 



C.l Optimum Length 



In order to use equation (4.24), intermediate values must first be determined 



j _ 0.00355 
“ (iV) 2/3 “ (0.7)2/3 



0.00450 



The temperature span parameter, r, cannot be calculated until T 5 is known, 
so an estimate based on previous iterations is used. If the calculated value 
of r is close to the value used, the computations are valid. 



{^/ r h)opt 



2 / r/N, t 

G \Bf{l/p c P c + l/p h P h ) 



( L/ r h)opt 



2 ( 0.365/0.00450 

15155 ( (0.287) (0.0078) ( 5 ^^+ n^J 



= 5283 
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Since r h is set by the geometry selected, L can be calculated 



L ovt = 5283(.0474/f)(.3048m//t) 
L opt = 19.08m 



C.2 Heat Exchanger Effectiveness and Pressure 
Drops 

From equation (4.16), 

Ntu, h = ( L/r n )N 3t 

N tu>h = (5283) (0.00450) = 23.79 
Applying the same assumption as was used in Chapter 3, 

N tU}h = Ntu,c — 2 Ntu 
N tu = 23.79/2 = 11.89 



Then, since this is still a balanced counterflow design, 



Nt 



€ = 



tu 



11.89 



= 0.922 



1 + Nt u 1 “I" 11-89 
For hot and cold side pressure drops, equation (4.17) is used. 
A P 

(—)* = f(L/r h )g*/(2pP) = f(L/r h )g* 



( AP ) 

\—p-)h 



= (0.0078) (5283)(0.05005) 2 = 0.010322 



(^)< 



= (0.0077) (5283) (0.00946) 2 = 0.00369 



A P 



= 0.10322 + 0.00369 = 0.1069 
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C.3 Entropy Generation 



Equations (4.14) and (4.15) will be used to determine entropy generation on 
both sides. 



Ns h = 



c *(&(¥). 



Ns h = 2£^ + (0.287)(0.10322) 



Ns h = 0.04497 



Ns e = 



01365 

23.79 



+ (0.287) (0.00369) 



Ns c = 0.0164 



Ns = Nsh + Ns c 



Ns = 0.04497 + 0.0164 = 0.0614 

C.4 Heat Exchanger Volume 

Since A/ r is calculated from minimum free flow area and other surface geo- 
metrical constraints, Af r from the corresponding case in Chapter 3 was used 
to calculate the heat exchanger volume. 

Vol = L(Afr) = 19.08m(0. 12876m 2 ) = 2.46m 3 

C.5 Cycle Efficiency Degrade 



Temperatures that reflect the actual effectiveness and pressure drop must 
first be determined. 

= Q-(l - k) = 5(1 - 0.1069) = 4.465 
Ps Pi 
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t 5 = 



P 4 \' , p‘( 2 T L ) 



t 4 



Ts = 



i (.OKt.4-1) 

(4.465) 1.4 

1300 



-l 



= 884.8°Jfc 



T 3 = e(T 5 -T 2 ) + r 2 



T 3 = 0.922(884.8 - 486.8) + 486.8 



T 3 = 853.8°Jfc 



= (r 4 - r 5 ) + { t \ - t 2 ) 

n t 4 -t 3 



v = 



(1300 - 884.8) + (300 - 486.8) 
1300 - 853.8 



rjoz ~ r) = 0.5760 - 0.5119 = 0.06412 
A check on the validity of the r value assumed shows 

■■ Wf -M 



f 1 884.8 i486. 

IV 486.8 “ V 88A 



0.367 



Therefore, assumed value of 0.365 was valid 



0.5119 
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Appendix D 

Sample Calculation of 
Entropy Generation for 
Conventional Designs 



Example uses case C-12. From previous example, 

= 0.34461 

= 0.0359 

= 0.00136 
= Ntu, c — 8 

Using equations (4.14) and (4.15), 




N, ih = + (0.287)(0.0359) = 0.05338 

Similarly, 

0 34461 

N SiC = — - + (0.287)(0. 00136) = 0.04347 

8 



APh 

Ph 

A Pc 

Pc 

Ntu,h 
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N s>h + N 3iC = 0.09685 



N, = 

Adding the first term of each side’s expression gives the AT contribution, 
and the sum of the second terms is the A P contribution to the entropy 
generation. 



N sAT = 0.08615 
N aA P = 0.01070 
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